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ABSTRACT 

We report the discovery of an active galactic nucleus (AGN) pair in the interacting galaxy 
system IRAS 20210+1121 at z = 0.056. An XMM-Newton observation reveals the presence of 
an obscured (N H ~ 5 x 10 23 cm -2 ), Seyfert-like (L 2 -iokeV = 4.7 x 10 42 erg s -1 ) nucleus in the 
northern galaxy, which lacks unambiguous optical AGN signatures. Our spectral analysis also 
provides strong evidence that the IR-luminous southern galaxy hosts a Type 2 quasar embedded 
in a bright starburst emission. In particular, the X-ray primary continuum from the nucleus 
appears totally depressed in the XMM-Newton band as expected in case of a Compton-Thick 
absorber, and only the emission produced by Compton scattering ("reflection") of the continuum 
from circumnuclear matter is seen. As such, IRAS 20210+1121 seems to provide an excellent 
opportunity to witness a key, early phase in the quasar evolution predicted by the theoretical 
models of quasar activation by galaxy collisions. 
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interactions favor both violent star formation as 
well as funneling of large amount of gas into the 
nuclear regio n to feed (and obscure) the accreting 



SMB H (e.g. JUrrutia et alj|200d ISchawinski et al 



2010t ) . T he importance o f merg e rs increases with 



redshift (jConselice et al] 120031: iLin et all 1.2008) 



and their fundamental role at the peak epoch of 
luminous AGN (i.e. quasar) and intensive star- 
formation activity at 1.5 z 3 is widely ac- 
cepted. 

Over the last few years, an increasing num- 
ber of interacting and disturbed molecular gas- 
rich galaxy systems showing both coeval power- 
ful starburst (SB) and quasar activity at high z 
have indeed been u nveiled (e.g., Carilli et al.ll2002l 
Dasvra et al.ll2008h . The counterparts in the local 



Universe to such luminous, high-z mergers are the 
(ultra)-luminous infrared (Ltr > 10 11 L Q ) galax- 
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ies, i.e. (U)LIRGs (|Sanders fc Mirabellll996h . In 
particular, these powerful objects should provide 
the opportunity for probing an inevitable outcome 
of the hierarchical merging process, i.e. the exis- 
tence of dust-enshrouded double/multiple SMBHs 
within the envelo pe of the host galaxy merger 
(jColpi et al.ll2008f ). Despite being widely pursued, 
direct observational evidence for AGN pairs in 
ULIRGs (as well as in all of the other types of 
galaxies) has been very limited so far. In the last 
few years X-ray observations with arcsec angular 
resolution have provided one of th e most efficient 
tools to disclose such systems. Komossa et al 



(|20031 ) discovered the first and unambiguous ex- 
ample of an active SMBH pair separated by d ~ 
1.4 kpc in the center of the ULIRG NGC 6240. 
Additional examples of dua l AGNs with a close 
sep aration were unvei led bv iBianchi et al. ( 20081 ) 
and lBallo etall (|2004l ) in the ULIRGs Mrk 463 (d 
~ 3.8 kpc) and Arp 299 (d ~ 4.6 k pc), respectively, 
on the basis of Chandra data. iGuainazzi et all 



(|2005af ) have reported the discovery of an X-ray 
bright AGN pair in ESO590-IG066, an early-phase 
(d ~ 10.5 kpc) merging system at z = 0.03. It 
is worth noting that in all these cases hard (>2 
keV) X-ray data have been crucial to detect ac- 
tivity from both SMBHs of the galaxy pair. In- 
deed, at least one pair member lacks optical/IR 
signatures of AGN activity, that suggests we are 
observing a non-standard AGN phase. 

Furthermore, a handful of kpc-scale dual AGN 
candidates have been recently uncovered in galaxy 
mergers by the detection of spatially-resolved, 
double-peaked emission line profile s with veloc 



ity of f sets of a few hun dreds km s 1 (jCivano et al 



20101 : iLiu et all |2010() . thus doubling the total 
number of bona fide active SMBH pairs collected 
so far and, in turn, opening interesting perspec- 
tives for further advances in this field of research. 

Here, we present the discovery of an AGN pair 
consisting of a Compton-thick (CT) Type 2 quasar 
and a heavily obscured Seyfert 2-like source in the 
interacting galaxy IRAS 20210+1121 (e.g., Sect. 
2). This discovery is based on the imaging and 
spectral analysis of XMM-Newton data described 
in Sect. 3. We discuss our results and conclude 
in Sect. 4. A cosmology with Ho = 70 km s _1 
Mpc -1 , r^A = 0.73 and fl\f= 0.27 is assumed 
throughout (|Spergel et al.ll2007l ). 



2. IRAS 20210+1121 

IRAS 20210+1121 is an interacti ng system of 
two galaxies separated by 12.2 arcsec (jArribas et al 
2004 iDavies et al.l l2002t ) . The larger component 
of the system, i.e. the southern galaxy (I20210S 
hereafter), shows a noticeable spiral arm structure, 
while the northern object (I20210N hereafter) is 
more spheroidal in shape. Furthermore, a bridge 
of emission connecting both galaxies is also visible 
in the optical band. However, narrow-band Ha 
imaging shows that I20210S has a centrally con- 
centrated, featur eless morphology, whil e I20210N 
is barely visible (jHeisler fc Vaderli 19951 ). I20210S 
is a LIRG (L IR = 7.8 x 10 11 L Q ) with a Seyfert 
2 nucleus at z = 0.056, and a [QUI] luminosity 
L[mm = 2.04 x 10 4 3 erg s" 1 (e.g., IPerez et all 



ill] 

19901 : IShu et al.ll2007l) 



I20210N has a very faint emission line spec- 
trum, but sufficient to derive the same z of the 
companion. An important feature of the op- 
tical spectrum of I20210N is an intensity ratio 
[Nil] A6584/Ha «3. Such a value is typical of both 
LINERs and Seyfert galaxies, and no firm conclu- 
sion on the presence of an AGN in this source can 
be drawn on the basis of these low signal-to-noise 
ratio data. Furthermore, a near-IR spectroscopic 
study of I20210N found a fea tureless near-IR con- 
tinuum (Bur ston et al. I l2nnih . 

I20210S was detected by BeppoSAX at a 2- 



10 keV flux level of 



2.9 



10 



-13 



erg 



cm s and only ^160 counts were collected 
in the 0.1—10 keV band (given the arcmin angu- 
lar resolution of BeppoSAX, any possible emission 
from I20210N cann ot be discerned). The analy- 
sis of these data by lUeno et ah (2000) revealed a 
very flat continuum slope T = 0.5 jl^'g and the re- 
markable presence of a strong (EWf = 1 -6ltf f 
keV) Fe Ka line that led these authors to sug- 
gest that I20210S may host a CT AGN. The ob- 
served F2-10 to [O III] flux ratio of < 0.1 in- 
ferred for I20210S also hints for a CT absorber 
scenario ( Guainazzi et al. I l2005bh . However, it is 
worth bearing in mind the large errors affecting 
the spectral parameters derived from this obser- 
vation. 
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3. Observations and Data Reduction 

We observed IRAS 20210+1121 with XMM- 
Newton on May 22, 2009 for about 75 ks (Obs. 
ID.: 0600690101). The observation was performed 
with the EPIC PN and MOS cameras operating 
in Full- Window mode and with the MEDIUM fil- 
ter applied. Data were reduced with SAS v9.0 
using standard procedures and the most updated 
calibration files. The event lists were filtered to ig- 
nore periods of high backg round flaring according 
to the method presented in lPiconcelli et al.l (|2004T ) 
based on the cumulative distribution function of 
background lightcurve count-rates and maximiza- 
tion of the signal-to-noise ratio. The PN source 
counts were extracted from a circular region of 
10 (I20210N) and 8 (I20210S) arcsec centered at 
(a 200 o = 20 h 23 m 25.04 s ; <5 20 oo = +11°31'47.7") 
and (a 200 o - 20' l 23 m 25.35 s ; S 20 qo = +11°31'30"), 
in order to avoid any cross-contamination between 
the two regions and include the maximum number 
of counts with E > 5 keV. For the two MOS cam- 
eras the extraction radius was of 9 arcsec for both 
sources. The background spectra were extracted 
from source-free, much larger circular regions on 
the same chip and close to the target. After this 
screening, the final net exposure times in case of 
I20210S(I20210N) were 61.2(60.1) and 69.7(70) ks 
for PN and MOS, respectively. Appropriate re- 
sponse and ancillary files for all the EPIC cameras 
were created using RMFGEN and ARFGEN tasks 
in the SAS, respectively. Spectra were rebinned so 
that each energy bin contains at least 20 counts to 
allow us to use the x 2 minimization technique in 
spectral fitting. 

In this Letter we present and discuss the PN 
spectral results only, since this detector has a bet- 
ter sensitivity over the broad 0.3-10 keV range 
compared to both MOS cameras (even if co- 
added together), and above 5 keV in particular. 
Nonetheless, we checked that consistent results 
were obtained including the MOS data in our 
analysis. 



4. Results 

An important result from the XMM-Newton 
observation is presented in Fig. 1, showing the X- 
ray images of IRAS 20210+1121 in the 0.5-2 and 
5-10 keV bands. In the soft band, only the X-ray 
emission centered on the [OIII]-luminous galaxy of 



the pair, i.e. I20210S, is clearly visible. Whereas 
the very hard X-ray image reveals the presence of 
two sources being almost comparable in intensity, 
and spatially coincident with the radio nucleus of 
the southern galaxy, or with the optical centroid 
of the northern galaxy (once an absolute astrom- 
etry uncertainty of 2 arcsec and the dispersion of 
photons due to the PSF are considered.]) , respec- 
tively. At their redshift, the projected separation 
between the intensity peaks of both sources in the 
5-10 keV image shown in Fig. 1 is d ~ 11 kpc. 
The simplest interpretation of the hard X-ray im- 
age is that I20210N may also host an obscured 
AGN, thus revealing, in turn, the presence of an 
AGN pair in this interacting system. Direct evi- 
dence to support this buried AGN pair hypothesis 
comes from the X-ray spectroscopy of I20210N. 

The spectral analysis of EPIC data of both 
sources was carried out using the XSPEC vl2 
software package. The Galactic column den- 
sity of Ng al = 9.7 3 x 1 20 cm" 2 derived from 
Dickey fc Lockmanl ( 199dh was adopted in all the 



fits. Henceforth, errors correspond to the 90% 
confidence level for one interesting parameter, i.e. 
Ax 2 = 2.71. 

We yielded a very good description of the 
spectrum of I20210N wit h a typical Comp ton- 
thin Scyfcrt 2 model (e.g.. iTurner et aL 1997 ). as 
shown in Fig. 2 (top panel). The primary X-ray 
continuum power law is absorbed by a column 
density of N H = (4.71} q) x 10 23 cm -2 and ex- 
hibits a slope of r = 2.0+0.2. The emission in 
the soft portion of the spectrum (the so-called 
soft excess component) is well fitted by an ad- 
ditional unabsorbed power law fixing its photon 
index to that of the absorbed power law, but with 
a different normalization (~ 3% of the primary 
continuum), plus three narrow Gaussian emission 
lines. The best-fit values for the energy of these 
lines are -0.82, -0.92, and -1.07 keV, which 
can be identified with Fe XVII 3d-2p, Ne IX Ka 
and Ne X Ka/Fe XXI 3d-2p transitions, respec- 
tively. Such a soft excess component has been 
detected in most of the obscured AGNs, and it 
is typically expla ined as emission from large-scale 
(—0.1-1 kpc; see Bianchi et al. ( 20061 )) photoion- 
ized gas, dominated by a wealth of strong emis- 



1 See http://xmm2.esac.esa.int/docs/documents/CAL-TN-0018.pdf 
for more detail. 
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sion lines from hydrogen- and helium-like ions of 
the most abundant m e tals, f r om carbon to sulfur 
dGuainazzi fc Bianchil 2007 : Kinkhabwala et all 
I2OO2J). Assuming this spectral model (x 2 /dof = 
15/20), we measured a 0.5-2 keV flux of 1.6 x 
10~ 14 erg cm" 2 s" 1 , and a 2-10 keV flux of 1.2 
x 10~ 13 erg cm~ 2 s _1 . After correcting for ab- 
sorption, this flux corresponds to a luminosity 
of L 2 -wkeV = 4.7 x 10 42 erg s _1 in the hard 
band. Such a value of the 2-10 keV luminosity 
falls well within the AGN luminosity range (e.g., 



5.2(6.6) x 10 41 erg s for the SB emission. Ac 



Maiolino et al.l 120031 ) , thus providing unambigu- 



ous evidence for the existence of an active SMBH 
at the center of I20210N. 

The XMM-Newton spectrum of I20210S is very 
complex as shown in Fig. 2 (bottom panel). This 
was expected on the basis of millimctcr/IR/optical 
data that have revealed the simultaneous pres- 
ence of star-forming and nuclear activity in this 
galaxy (Horellou et al. 19951 : Burston et al. 2001 : 
Perez et al.l Il990h . In particular, from the 1.4 



GHz radio(far-IR) luminosity of I20210S, a star- 
formation rate SFR ~ 120(75) M0 yr^ 1 can be 
es timated according t o the relationship reported 



Ranalli et al 



(2003). We have indeed found an 
excellent description of the XMM-Newton data 
assuming a composite SB + AGN emission model 
( x 2 /dof = 0.96(65)). The soft X-ray SB emission 
has been fitted by the superposition of two thermal 
emission components (MEKAL model in XSPEC) 
with solar metallicity and a temperature kT = 
0.58±0.08 and kT = 1.25±£^ keV, respectively, 
in agreement with typical values of temperature 
mea sured in other we ll-known star-forming galax- 
ies (jPtak et al.l 1 1999T) . The hard X-ray emission 



has been described assuming an absorbed cutoff- 
power law model in the form E~ T exp~ ,w / kT with 
photon index T = 1.1±0.2 and cutoff energy fixed 
to 10 keV. Such a spectral shape is expected in 
case of a contribution from flat-spectrum bright 
Low-Mass X-ray binaries and, mostly, from High- 
Ma ss X-ray binaries (HMXBs), as pointed out 
by |P ersic fc Raphaelil (|2002h . Another striking 
characteristic of the emission from HMXBs is a 
strong Fe XXV emission line at 6.7 keV, that 
is indeed observed in the spectrum (see Fig. 2, 
bottom panel) with a poorly-constrained value 
of equivalent width EW ~ 500 eV. We derived 
a soft(hard) X-ray luminosity of -Zjo.5-2fceV = 



cording to iRanalli et al.l ([2003) , these values im- 
ply a star-formation rate SFR ~ 110—130 M Q 
yr , which is consistent with the SFR values 
derived both from the radio and far-IR luminos- 
ity reported above. The goodness of this match 
therefore lends further support to the idea that 
the hard X-ray power-law component originates 
from the population of HMXBs expected to be 
present in the SB regions of 12021 OS. 

Our best-fit model to the XMM-Newton data 
includes a component due to X-ray reflection from 
cold circumnuclear matter with Nh £ 1.6 x 10 24 
cm~ 2 (i.e. CT), and an Fe Ka emission line to 
account for the reprocessed AGN emission visible 
in the 0.3-10 keV band, the X-ray primary con- 
tinuum emissi on being completely b locked in this 
energy range (jGhisellini et al.lll994l ). 



The energy centroid of the line is at 6.351q'q4 
keV and the EW measured with respect to the 
reflection continuum is 900±400 eV. These val- 
ues unambiguously indicate an origin from re- 
flection in cold circumnuclear CT material. The 
X-ray primary continuum from the AGN is to- 
tally depressed in the EPIC band, as expected in 
case of a CT absorber. We derive an AGN flux 
(F2-10 = 7.7 x 10~ 14 erg cm~ 2 s^" 1 ), accounting 
for -47% of the total (AGN+SB) 2-10 keV flux. 
The observed L2-iofceV of the AGN component is 



5.3 x 10 41 erg 



according to the reflection- 



dominated/CT scenario it should be a t most 1- 
2% of the de-absorbed lum inosity (e.g.. IComastri 
120041 : iLevenson et al.ll2006l ). suggesting an intrin- 
sic I/2-iofceV > 0.5-1 x 10 44 erg s -1 , which is con- 
sistent with the expectation based on the -Zjjoiii]- 
This result matches well with the hypothesis of a 
CT absorbing screen and, hence, the presence of a 
quasar 2 (with Lx >10 44 erg s _1 ) at the heart of 
I20210S. 

We also tried a model assuming a transmission 
scenario for the AGN emission below 10 keV. We 
fixed the photon index of the continuum power law 
to the canonical value of T = 1.8 due to the lim- 
ited statistics. This fit is statistically as good as 
the reflection-dominated fit discussed above, with 
a resulting N H = (3.2± 2 ;°) x 10 23 cm" 2 . However, 
this Nh implies an L 2 -iofceV = 1.2 x 10 42 erg s _1 , 
which is two orders of magnitude lower than ex- 
pected on the basis of the £[om], an( A an EW of 
the Fe Ka line against the absorbed continuum of 
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~ 250 e V (|Ghisellini et al.l Il994t iGuainazzi et al 
2005bl ). while we measured an EW = 620±260 
eV. These two considerations tend to disfavor the 
transmission scenario and lead us to assume the 
presence of a CT screen along our line of sight to 
the nucleus of I20210S as the most likely interpre- 
tation of the XMM-Newton data. 

5. Conclusions 



The XMM-Newton observation of IRAS 20210+1121 
presented here has unveiled the existence of an ob- 
scured AGN pair placed at a projected distance 
of d ~ 11 kpc in this interacting galaxy system. 
In particular, we have discovered a Seyfert 2-like 
AGN in the nucleus of I20210N, for which nei- 
ther optical nor near-IR spectroscopic observa- 
tions have provided unambiguous evidence for the 
existence of an active SMBH at its center. Fur- 
thermore, the results of our spectral analysis have 
provided evidence that the southern member of 
the pair, the LIRG I20210S, optically classified as 
a Seyfert 2 galaxy, likely hosts a powerful AGN 
hidden behind a CT absorber. The AGN is embed- 
ded in a strong SB emission accounting for ~50% 
of the 2.0-10 keV flux measured for I20210S. 

IRAS20210+1121 is therefore a rarely-observed 
example of CT quasar 2 plus optically 'elusive' 
AGN pair observed during the initial stage of 
the interaction between their host galaxies, that 
are still easily identifia ble, but also show a well- 
developed tidal bridge (jArribas et al] 120041 ) . As 
such, this system seems to provide an excellent 
opportunity to witness a merger-driven phase of 
quasar fueling predicted by most of the evolution- 
ary models based on the co-evolution of SMBHs 
and their host galaxies. 

The mismatch between the optical/near- IR and 
the hard X-ray appearances of the nuclear spec- 
trum of I20210N can be explained in terms of a 
completely blocked line of sight to the nuclear re- 
gion, so th at the narrow l i ne reg ion (NLR) is also 
obscured (|Maiolino et al.l 120031 ). The geometri- 



cal properties of the absorber should be thereby 
different from those assumed in typical Seyfert 2 
galaxies, i.e. a pc-scale torus-like shape. The ab- 
sorber in I20210N may be characterized by a much 
more extended distribution (over a few hundreds 
of pc) in a way that obscures the NLR. Alterna- 
tively, it could be spherically symmetric blocking 



the flux of ionizing UV photons responsible for the 
line emission in the NLR. The AGN in I20210N 
is revealed through hard X-ray observations only; 
this feature is shared by most AGN pairs dis- 
cover ed lately thanks to hard X-ray observations 
(e.g., Komossa et al. 20031 IGuainazzi et al.l 2005a ; 



Bianchi et al.l 120081 ) . An intriguing explanation 
for this behavior can be given in terms of a 
dust-enshrouded circumnuclear environment due 
to merger-induced processes favoring gas concen- 
tration in the galaxy center. For instance, there 
may be an extra-torus optical/X-ray absorber ly- 
ing far from the nucleus, and outside the NLR, 
being likely associated with prominent dust lanes 
in the disturbed host galaxy. 

The discovery of a CT quasar 2 in I20210S is in 
itself very important given the paucity of low— z 
members of this peculiar class of AGN detected 
so far, which are considered a key ingredient in 
the synthesis models of the Cosmic X-ray Back- 
ground. This can be ascribed to their low surface 
density and their absorption-induced faintness at 
the wavelengths where classical large-area surveys 
have been performed (i.e. optical, near-IR, UV, 
X-rays), that make the luminous CT AGN pop- 
ulation extremely difficult to observe. Selection 
criteria based on m id-IR vs. optical colors (e.g., 
Lanzuisi et al.l 120091 and references therein) have 
been proven to be efficient in discovering a large 
number of heavily obscured quasar candidates at z 
> 1. Unfortunately, most of these sources are de- 
tected at very faint flux levels (F 2 _io <C 10~ 14 erg 
cm~ 2 s _1 ), making an appropriate X-ray spectral 
follow-up extremely time-consuming. 

Further support to the CT quasar 2 nature for 
the AGN in 12021 OS is provided by the compar- 
ison of the expected value of the X-ray bolomet- 
ric correction rx,bol = ^2-iofcev/Lboi = 0.043 x 
(Lboi/10 45 ) -0 ' 357 (assuming a bolometric luminos- 
ity of L bo i(« L T r ) w 3 x 10 45 er g s" 1 ) for a typ- 
ical quasar from IPiconcelli et al. (20071), and the 
values of rx,boi calculated using the hard X-ray 
luminosity inferred for the CT and transmission 
scenario, respectively. In fact, the expected value 
of rx,boi ~ 0-03 is consistent with that derived as- 
suming a L2-wkeV > 5 x 10 erg s" 1 (estimated 
from the £[oin] luminosity and a CT absorber), 
i.e. rx,boi > 0.02, but it is much higher than rx,boi 
= 0.0004 derived for a L 2 -iofc e v = 1.2 x 10 42 erg 
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Detection of objects such as I20210S is impor- 
tant because they provide useful templates to ex- 
plore the multiwavelength properties of the ob- 
scured accretion phenomenon without any lumi- 
nosity bias. The observed properties of I20210S 
can be interpreted in the framework of an evolu- 
tionary merger-driven scenario according to which 
a peculiar dust-cocooned, early stage in the life 
cycle of quasars is linked to a period of intense 
star-forming activi t y in the interacting h ost galaxy 
(|Silk fe Reesl 119981 iTreister et al.l l2010h . An eas- 
ily observable outcome of this scenario is indeed an 
enhancement of the IR luminosity, with the system 
undergoing a (U)LIRG phase po wered both by 
the SB and the AGN. Furthermore. iHorellou et"al 



(|1995I ) measured a molecular hydrogen mass Mh% 
of 4.1 x 10 9 M Q for I20210S that implies an 
Lfir/Mh2 ratio > 100 L^M^ 1 , i.e. a value typi- 



cal for gas-rich mergers (panders et al.lll991l ). Ac- 
cording to model predictions, the SMBH at the 
center of I20210S should be accreting close to the 
Eddington rate, in agreement with the quasar-like 
values of £[om] ) measured(estimated) 

for the AGN. 

The interacting system IRAS 20210+1121, 
therefore, surely deserves deeper investigations in 
the future in order to examine the possible pres- 
ence of any structures associated with the merg- 
ing process (i.e. outflows, inflows, obscured SB 
regions) which could not be revealed by the obser- 
vational data available so far, but potentially very 
useful for our understanding of quasar evolution 
and AGN/SB triggering mechanisms. 
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Fig. 1— XMM-Newton EPIC PN Gaussian- 
smoothed (<r — 2 arcsec) image of the interact- 
ing galaxy system IRAS 20210+1121 in the 0.5-2 
keV (left panel) and 5-10 keV (right panel) energy 
range. 
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Fig. 2. — Top: Best-fit data and folded model 
(thick line), plus residuals, of the PN spectrum of 
I20210N. Bottom: The PN spectrum of I20210S 
with the "composite" CT AGN+SB best-fit model 
plotted as a thick line. The individual model com- 
ponents are also shown: (a) pure reflection con- 
tinuum resulting from a power law illumination 
of cold material; (b)+(c) thermal plasma compo- 
nents; (d) cutoff power-law associated with the X- 
ray binaries; and (e) two narrow Gaussian lines at 
6.4 and 6.7 keV, respectively. See text for details. 
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